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ABSTRACT 


This report contains four parts of the Reynolds-stress closure modeling: 
1) Improvement of the k and c equations, 2) development of the third-moment 
transport equation, 3) formulation of the diffusion coefficient of the momen- 
tum equation by using the algebraic-stress model --df turbulence, and 4) the 
application of the Reynolds-stress model to a heat exchanger problem. It has 
been demonstrated that the third-moment transport model Improved the predic- 
tion of the triple-velocity products in the recirculating and reattaching flow 
regions In comparison with the existing algebraic models for the triple- 
velocity products. Then, optimum values for empirical coefficients are ob- 
tained for the prediction of the backward-facing step flows. Furthermore, a 
functional expression Is derived for the coefficient of the momentum diffusion 
by employing the algebraic-stress model. Finally, the second-moment closure 
is applied to a heat transfer problem. The computations for the flow in a 

l 

corrugated-wall channeJ show that the second-moment closure Improves the pre- 
diction of the heat transfer rates by 30X over the k - c model. In this 
way, the application of the Reynolds stress model to heat transfer problems is 
demonstrated to be promising 



PREFACE 

This report outlines the progress of the research conducted from July 
1905 to January 1986. During this period, the third-foojccnt closure model was 
studied as well as refining the second-moment closure. The development of 
transport equations for the third-moments of turbulence velocity fluctuation 
has been pursued since It was observed. In the last research project period, 
that none of the existing algebraic third-moment models can predict the levels 
of the triple velocity products correctly for the backward-facing step flows. 
The predictions by the algebraic models were particularly poor In the recircu- 
lating flow regions. 

In this report It is shown that the currently developed transport equa- 
tions of the third-moments Improve the prediction of the triple-velocity pro- 
ducts of turbulence fluctuations In comparison with the algebraic third-moment 
models. 

In carrying out this project the authors are grateful for the assistance 
of Messrs. A. Bagherlee, R. Smith, T. Niess. and V. Kodali who performed a 
large amount of programming and computations. 
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NOMENCLATURE 

a 

channel step height 

4 a 

a u 

anisotropy 

I 

b 

channel width 


C l* C 2 ,C p* ' 

C k’ C p’V 

C wl ,C w2 ,C c* 

i constants used In turbulence model 


C c1 ,C c2 ,C <,r 

C «.2 ,C $3* 


4 

C 1J 

coefficients for near-wall Reynolds stresses 

. 

D 

diameter of the disk 

1 ’ i 

°U 

diffusion rate of the Reynolds stresses 


°r 

diameter of the tube 

■ | 

f 

average skin friction factor 

i 

f w 

function for wall correction 

\ . 

6 

gene *at Ion rate of turbulence kinetic energy, k 

’ I 


generation rate of Reynolds stresses 

! | 

h 

local heat transfer coefficient 

* 1 

H 

step height 


H U 

secondary generation rate of Reynolds stresses 

* 

k 

turbulence kinetic energy (=u^) 


L fin 

length of fin 


Nu 

Nusselt number based on 2b 


Hu 

averaged Nusselt number 

< 

P 

pressure fluctuation 

1 

! 

P 

mean pressure 


<*w 

wall heat flux 

1 

r 

radial coordinate 

l 

l 

i 

\ 

i 

• - ’ r 

\ L 





x, y. z 


Reynolds number based on 2b 
temperature 
bulk temperature 
wall temperature 

fluctuating velocity In x direction 

mean velocity 

Inlet stream velocity 

averaged mean velocity 

fluctuating velocity In y or r direction 

mean velocity In y or r direction 

fluctuating velocity In 2 or 0 direction 

Cartesian coordinates 


* tt 3 


dimensionless distance from wall to the first numerical 
1/2 

node-point, (- k y/v) 
constants used in £q. (19) 
thermal dlffuslvlty 
Kronecker delta 
energy dissipation rate 
coordinate in azimuth direction 


dynamic viscosity 
kinematic viscosity 
density 


Prandtl number 


pressure-strain correlation 


Subscripts 


tensor notations 
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1 . INTRODUCTION 

In computations of turbulent separated shear flows. It has become more 
common to use the second-order closure of turbulence. While modem computa- 
tional techniques enable us to successfully predict most of the simple turbu- 
lent shear flows. It Is still difficult to accurately predict the separated 
and reattaching shear flow particularly when It Is accompanied by a flow 
recirculation. This Is partly because the computation should be performed 
with an elliptic approach which requires large amounts of computer time and 

storage, and partly because a turbulence model that can predict such complex , 

* 

turbulent flows (Including flow reattachment, flow recirculation and a t 

recovering boundary layer) has not yet been completely developed for universal 
usage. ' 

A large number of experimental studies on this subject have been reported ! 

In the last two decades. The Reynolds stresses In the reattachment zone were j 

12 3 I 

obtained by Etheridge and Kemp , Kim et al. , Smyth , Eaton and J 

4 

Johnston , etc. Commonly It has been shown that the turbulence energy level 
reaches a peak value approximately one step height upstream of the reattach- f 

ment point, then decays rapidly In the streamwlse direction toward the wall, 
although It decays relatively slowly along the wall in the wall vicinity 
region. This feature Is In contrast with free shear flows which are not con- 
strained by the solid wall. It Is also noted that the recirculating zone 
affects the reattaching shear layer resulting in higher turbulence energy 
levels near the reattachment zone. 

Although the measurements of nigher-order turbulence moments are scarce 
In literature, a few third-moment data are available. The triple-velocity 

5 

products were measured by Chandrsuda and Bradshaw by using a hot-wire 
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anemometer for a channel with a larger channel width of 3.5 step heights. 
Driver and Seegralller 6 used a channel with a width of 9 step heights wherein 
the boundary layer thickness Just above the step was about 1.5 step heights 
having an Invlscld core region outside the separating layer. They measured 
the secondand third-moments of turbulence velocity fluctuations with a two- 
component IDV by varying the pressure gradient In the Invlscld core region. 

Through this NASA project, the authors have Investigated the features of 
the second- and third-moments In the reattaching and recirculating flow re- 
gions. First, the existing model of the Reynolds-stress transport equations 
was modified and adjusted for the elliptic flows. 7 The computations were 
made for several different cases of flow and geometry, and the results were 

g 

compared wltn the k - c model and the algebraic-stress model of Rodl. 

Although the model used by Amano and Goel 7 was the hybrid model of the 
Bousslnesq viscosity and the Reynolds-stress equations, the prediction of the 
normal and shear stresses was improved In the reattaching and the recovering 
region behind a s*ep. 

Meanwhile, It was observed that the triple-velocity products decay 

rapidly In the reattaching region which creates steep changes In the diffusion 

9-10 

rates of the Reynolds stresses. In the last two reports, most of the 
existing third-moment models were tested to determine whether they could be 
employed for the prediction of such reattaching shear layers or not. As a 
result. It was recommended that a transport model for the third-moment should 
be evolved since none of the algebraic models for the third-moments gives 
sufficiently high levels in a complex turbulent flow although the models are 
satisfactory for parabolic thin shear flows. 


In this report, the work performed Is sunsaarlzed In the following steps: 

1. Incorporation of the non-lsotroplc treatment for the diffusion rates 
of the k and c equations. 

2. Determination of the empirical coefficients of the above mentioned 
diffusion rates. 

3. Development of the transport equations for the third-moments of 
turbulence. 

4. Testing the third-moment model for the reattaching shear flows In 
order to determine the coefficients which appear In the third-moment 
transport equations. 

5. Formulation of the correlation function used in the diffusion rates 
of the momentum equations. 

6. Application of the second-moment closure to a heat transfer problem. 
In the above procedure, step 5 was achieved by using the algebraic-stress 
model because of its simple formulation of the Reynolds stresses. In order to 
eliminate the complexity due to the wall flows, this test was limited to only 
wakes which were created behind steps. In step 6, a periodically corrugated 
wall channel was considered. This Is the basic study of a compact heat ex- 
changer. 


2. THE SEC0ND41CHENT CLOSURE KODEl 
2.1 Prellmlna.v Remarks 

Both the turbulence energy, k, and the energy dissipation rate, *, are 

3/2 

always used to define the time scale, k/e, and the length scale, k /c. 
Therefore, If the k- and c-equatlons are not functioning with the stress 
equations In the same order of treatment, the modeling Is not considered to be 
complete, and thus the overall results will not be Improved. To deal with the 
higher order forms of these equations, the original (unclosed) forms need to 
be employed to define the turbulence energy and the energy dissipation rate 
and then some empirical correlations have to be reevaluated for the flows 
considered In the present study. 

The generation and dissipation rates of the k-equatlon do not require any 
empirical coefficient. On the other hand, the generation and dissipation 
rates of the c-equatlon have empirical coefficients which have been used by 
a large number of researchers for more than a decade Therefore, the key 
terms which may improve the model are the diffusion rates of the k- and 
c-equatlons that need revising for the second-moment closure. 

In the second-moment closure, the conventional Bousslnesq viscosity (dlf- 
fuslvlty) model should be replaced by the second-moment expression. In this 
process, the empirical coefficient for the diffusion rates of the k- and 
c-equatlons must be determined. However, these coefficients have not yet 
been evaluated for universal usagi*. The determination of the values for these 
diffusion coefficients Is attempted In this section for the better prediction 
of the reattaching and recirculating flows. 
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2.2 Improvement of k. c and the Reyno^ds-Stress Equations 

The transport equation which describes the Reynolds-stress variation in 
an elliptic flow field is derived as 


, !Vi 

J i ax. 


3U 


1 


-< Yi + u < u 

(i) 

au.u. 


au 
i u i ax 


„ au. au. 

> <sr + sr 1 ) 

l p ax J 2X 1 
01) 


au, au. 


, OUtUi VUj ou. 

+ Vi + Vj> - VjY • 2w 

(111) (1v) 


(i) 


from which the turbulence kinetic energy Is also derived. The definition of 
the turbulence kinetic energy is 


k = u*/2 (2) 

2 2 

There are two ways to obtain k; one is to solve Eq. (1) for u , v , and 
2 

w and substitute the results into Eq. (2), and the other way is to formu- 
late the transport equation for k itself. When the problem is In a two- 
dimensional plane coordinate system, the number of equations which need to be 
solved is the same either way. Here the transport equation for k is converted 
Into the second-moment closure rather than solving all three normal components 
of the Reynolds stress and summing them to obtain k. This is because the 
k-equatlon is widely used as an Indicator of turbulence levels. Thus, the 
k-equatlon cannot be omitted when one is solving the two-equation model or the 
algebraic-stress model. 

Upon contracting Eq. (1) by setting j equal to i and dividing by 2, it 


yields the k-equatlon as follows: 


tttt gU i _ a , aic 
u l u a ax a ax t ax t 


rr . i!v _ 2i 

P *■ 2 ^ v 3x. 3X. 

t & 

(11) (111) 


Terms (1) and (111) on the right-hand side of Eq. (3) represent the generation 
and the dissipation rates, respectively. The dissipation can be determined by 
solving the c-equation. Term (11), the diffusion rate, may be replaced by 
the following form. 


(11) - 


a *t ^ c Ul ° l Sx i^ 


where it is assumed that the second term In the parenthesis is negligibly 
small in comparison with the third one In (11) In a large-scale turbulence. 
This assumption Is supported by Irwin's measurements of self-preserving wall 
jets. 11 The first term in (11) Is kept in the computation in order to 
account for the viscous effect near the wall. 

To obtain the Reynolds stresses. Eq. (1) was closed as follows: 

3 V7 3U. 3U< 

U t ax 4 a - (u j u ! 3^ + Vi + *ij + Dij - f V 

where the pressure-strain correlation. is divided Into three parts: 
the Rotta term 12 , the rapid term 13 and the viscous term. 14 These are 
combined and given as follows: 


>1J " " C *1 k (U i U j " 3 \i - C 62 (G 


'♦2 ' ij 


3 V> 


+ C «3 {H 1j “ 3 4 1j 6 > 


- Q - 


a l 



where 


au. 


G “ - u,u 


axj 


au 


J U 


au. 


(u 1 u kS^ + Vk a^> 


au.. 


au,. 


H 1J " * (u i°k SJ + u j u k 


(7) 


C *1 - 1.5 - 0.125 f w 


C *2 - °' 4 ‘ °- 015 f w 


C *3 - 0.015 f w 


and where f y is a function that controls the strength of the wall 
correction. Here it is taken to be proportional to the length scale. 


k 3/2 /c. 


The diffusion rate, D^j, is conventionally approximated as 


a l. au.u, 

D U • it < c d ) v. 

«• m 


( 8 ) 


for the computation of u^Uj which is a rather rough approximation because 
the diffusion rate contains the triple-velocity products. A more extensive 
model for the triple-velocity products u^u^, is given in the next 
chapter. 

The energy dissipation rate :qi:ation is given as 



I 


J 

1 ax. 


, su* 

- k (C-1 T~ * c 


'cl i £ a*. 

I 

0) 


3u^ au 1 

« * 

tit) 


ax. 


(w u 


a “i?X 2 sM !!i n_. 

- a< j "j ” «, 


t sx a ax ft 


(111) 


(1v) 


(V) 


( 9 ) 


where C cl «* 1.45 and C f2 - 1.92 which have been used to date by a num- 
ber of researchers as was discussed In section 2.1. The generation (1) and 
the destruction (11) terms are evaluated using a direct approach. The diffu- 
sion rate Is modified In a similar fora as Eq. ( 4 ) : 


(ill) + (1v) + (v) 




) 


( 10 ) 


The coefficients C fe and C ( were recommended to be 0.31 and 0.15, respec- 
tively. by Pope and Whitelaw. 15 since these values are appropriate only for 
simple free shear flows, more extensive tests for these coefficient-, are made 
for the computations of reattaching shear flows. The optimum values for these 
coefficients are recommended In the following subsection. 


2-3 Computed Results and Discussion 

The flow region considered is given In Fig. 1. At the inlet of the flow 
field the prescribed values are given to all variables; they are be .Ically 
taken from the experimental data with which the results are compared. Along 
the top portion and the outflow section of the computational domain, the 
contlnuatlve boundary condition Is applied. 
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At the wall boundaries, the “wall law" Is used to specify Rean veloci- 
ties, turbulence kinetic energy, and the Reynolds stresses. On determining 
proportional constants between k and In the wall-vicinity region, we 
obtained the following relations for u^: 




<i-*14> 


Y 

p dx 


OD 


where x and y are stroamwlse and transverse coordinates, respectively. The 

coefficients C,. used here ore taken from' Ref. 9. 

IJ 

C n - 1.21, C 22 - 0.24, C 12 « 0.24. (12) 

The dissipation rates along the wall are obtained under the “local equilibrium 
condition*; that is, the length scale near the wall Is proportional to the 
distance from the wall. 

The solution domain consists of a 60Hx5H area with 52x52 grid points. 

The grid expands at rates of 3X and 2X In x and y directions, respectively. 

This system was selected after a number of grid tests ‘o produce the optimum 
grid-independent condition for the computation of mean velocities and the 
second-moments of the turbulent velocity fluctuations. 

Figures 2 through 6 show the computed results of mean velocity and the 
Reynolds stresses at several streasiwlse locations behind the step. By using 
the present model for Re^ - 32,000, the dimensionless distance, y + , from 
the solid wall to the first numeric; 1 node-point varies from 16 to 37. Here 
Re^ denotes the Reynolds number based on the step height H. The typical CPU 
time consumed for one computation took approximately 30 minutes on a UNIVAC 
1100 with about 300 Iterations. It Is also found that the CPU time depends 
slightly on the values of the coefficients C^, and C £ In Eqs. (4) and 
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'I 

•i 

(10). Maaoly, the computations with •» 0.31 and ■ 0.15 took about 
16% more CPU time and 50 additional Iterations than In the cases of - 0.1 

and C • 0.3. 

c 

As shown In Figs. 2-6, the profiles of the Reynolds stresses computed 
with the coefficients C^ - 0.31 and - 0.15, which are presently 
recommended for free shear layers, are about 20% to 30% lower than the experi- 
mental data of Driver and Seegmlller 6 for the normal stresses, and 14% lower 
for the shear stresses. As a result, the levels of the mean velocities In the 
shear layer region are 20% higher than the experimental data due to a reduc- 
tion In turbulent diffusion rates. 

In order to Improve the prediction of these turbulence stresses in the 
reattaching shear layer, a number of parametric tests have been performed for 
different values of the diffusion coefficients. It was discerned that the 
smaller values for and larger values for C £ gave better results in 
comparison with measured data. This Is because the larger values of 
cause the diffusion of turbulence energy to Increase, thus resulting in lower 
levels of the Reynolds stresses. Similarly, the smaller values of C e , In 
turn, increase the dissipation rate of the turbulence energy which results In 
a reduction In the turbulence energy. 

i 

In the present flow field an extra effect Is created by the recirculating 
flow which occurs below the reattaching shear layer. The recirculating flow 
at Ihe corner of the step and the bottom wall enhances the level of turbulence 
energy in the separated flow resulting In much higher energy at the reattach- 
ing region, and this energy is transported downstream. 

After performing parametric tests It was found that the combination of 
C k ■ 0.1 and - 0.3 gave the best results for both the mean velocity 
and the Reynolds stresses. It Is also observed In Figs. 2-6 that even the 
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Bousslnesq sodel* gives better results than the Reynolds-stress model with 
C fc • 0.31 and » 0.15. The Improvement fade by the Reynolds-stress 
model with the new coefficients C k - 0.1 and C c » 0.3 Is primarily at- 
tributed to the Inclusion of the second-moments In the diffusion terras which 
account for the nan-isotropic effects. Table 1 compares three models 
discussed above. 



A 


1 ! 


TABLE 1. Diffusion terms for k and e equations 



Ck 

C c 

c„ 

Agreement With 
Experiments 

Reynolds-stress model 
Pope and Whltclaw 
Present Study 

0.31 

0.1 

0.15 

0.3 

Not Used 
Not Used 

Bad 

Good 

Bousslnesq Model 

Not Used 

Not Used 

0.09 

Fair 


The turbulence energy balance Is examined by using both the presently 
determined Reynolds-stress model and the Bousslnesq model (cases 2 a,;d 3 In 
Table 1, respectively) and Is shown in Figs. 7-11. The variation of the 
convection terra (Fig. 7) shows that both models give equally reasonable levels 
at several different locations. For the distributions of the diffusion terra, 
however, the computed results do not agree well with the experimental data in 
the recirculating region, as shown at x/H - 3.9 In Fig. B. This Is the reason 
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« 

I 

j 


1 * • 
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*The Bousslnesq model has the diffusion terra In the following form: 

\ • 3 ^ <°-°’ £ 

where $ represents k or c, and o stands for Prandtl number for k or c. 
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why the triple-velocity products that represent the diffusion of the Reynolds 
stresses are examined by using the transport equations for UjUjUj, In 
this project. Botn the computed distributions of the production (Fig. 9) and 
the dissipation (Fig. 10) are relatively In accordance with the experimental 
data although both models always slightly overpredict these quantities. 

The overall balance of these terms Is compared in Fig. 11. Here the 
computations are made by using the Reynolds-stress model with the presently 
recommended diffusion coefficients. It Is observed that both the production 
and dissipation rates predominate In the shear layer near the step, but these 
levels decay quickly toward downstream. On the contrary, the diffusion and 
the convection rates stay almost at constant levels In the streamwise direc- 
tion while they vary rapidly In the transverse direction. 

2.4 Summarizing Remarks 

In this section. It has been shown that both and C e are con- 
siderably Influential in the results of turbulence qu intitles and also In the 
numerical convergence history of the Iteration procedure. This Indicates that 
the diffusion process of the turbulence energy In the reattaching shear layer 
Is significant which Is in contrast with the boundary layer flows. Thus, It 
Is concluded that the determination of the diffusion coefficients should be 
done carefully. Probably, a more elaborate functional expression will suffice 
the requirement rather than employing the constants. 

Secondly, it became clear that the third-moment closure should be Intro- 
duced for better understanding of Ihe diffusion process In the flow recircula- 
tion region. This is carried out in the next chapter. 


3. THE THIRC-HOfiEHT CLOSURE KODEL 




3.1 Preliminary Remarks 

In closing the Reynolds-stress equation (1), the terms of the pressure- 

strain (11) and the diffusion (111) have been approximated by Eqs. (6) and 

(8), respectively. The question arises whether these approximations are 

universally accepted or not. With respect to the pressure-strain correlation, 

several models have been examined In the latest report, 1 ® and the tests for 

13 1 & 

the best model revealed that three models by Naot et al., Naot et al., 

14 

and Launder et al. produced essentially the same results. In this 
chapter, attention Is primarily focussed on the diffusion term (term (111) In 
Eq. (D). 

3.2 Existing Models and the Transport Model for the Third-Homents 

Terra (111) In Eq. (1) represents the diffusion of the Reynolds stresses 
due to molecular viscosity, pressure fluctuations, and turbulence velocity 
fluctuations. It has been shown that the diffusion due to molecular viscosity 

and pressure fluctuations is negligibly small at large turbulence Reynolds 

* 

numbers. Thus, the triple-velocity products are the most dominant diffu- 
sive agency for the Reynolds stresses. 

It Is known that the higher moments of turbulence fluctuating velocity 
vary rapidly In the reattaching shear layer towards a solid wall.^ This 
affects the accuracy of the computation In the diffusive term of the Reynolds 
stresses because It contains the third-moment velocity fluctuations. It was 
discussed in the last report 10 that the triple-velocity products need to be 


♦'These terms appear to be small In most circumstances since turbulence energy 
budgets balance, to within experimental error, even when these are neglected. 
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evaluated In a transport equation model rather than In on algebraic fom be- 
cause the convection and generation of the triple-velocity products have to be 
taken suitably into account according to the change in the mean strain rates. 
For this reason a full transport equation model of the triple-velocity pro- 
ducts Is developed and used for the prediction In the reattaching shear flows. 

The complete equation is given as 


u i <"t “jV - - <Vj u i sj + u j Vi *7 * vi"» «;> 


a “k“l 


a Vi 


!!!&> 


* <u i“j * “A ♦ Vt > 


- 3^ ‘VAV - 1 <Vj S 7 + “A S 7 ♦ Vi Sj> 


<iii) 


where terms (1) and (11) represent the generation due to the mean strain rate 
and the generation due to turbulence stresses, respectively. Terms (ill) and 
(1 v) In Eq. (13) both represent the diffusion rate of the triple-velocity 
products. Here term (ill) Is correlated by assuming Gaussian form as follows: 


cun - — < Vj . Vl ♦ Vk • Yi * Vj ‘ Vi> 


and term (iv) Is given as: 


(iv) -- C p jru lUj u k 
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Then the final form of Eq. (13) becomes 


J i (U 1 Yk> " - < u i U A 5^ + u j V* ixj + Vi u a JT) 


a au i u i u k 
ax. v ax. ' 


“ (u t u, 


au 1 u k 


t“* "«r + u -“‘ -=r * v. -sr» -s f w 


In the discussion section the results of the parametric tests are shown for 
optimum values for the coefficient C p . The term with the laminar viscosity 
v Is added to account for the viscous effect. 

The above mentioned model Is also compared with four existing algebraic 
models. These are given as follows: 

(1) The model of Daly and Harlow 17 


u 1 u ju k * 


k aujU. 

- °- 25 r vt -sr 1 

t 


(11) The model of Hanjallc and Launder 18 


k au.u. au.u, 

U 1 U j u k “ - 0* 11 c I Vj 17~ + Vi + Vk *aB 


(111) The model of Shir 19 


2 au.u 


Vj u k - - °- 04 r » 
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. j*JL 


1 


( 1 v) Ths model of Corraack et el. 


20 


u 1 u j u k “ 4 C {2< *1 (4 ij 4 kt + 4 ik 4 ja * 4 kj 4 U } 


ak 

ax. 


+ °2 (a 1k,J + a 1J„k + a kj,1 )} 


+ 2 < l2 “3 (i 1k a Jl * ‘tj a kl + *jk s U> W. 


(19) 


3.3 Computed Results and Discussion 

In order to establish a reliable model for the prediction of the thlrd- 
f.oments in the reattaching shear layer, the values of the mean velocities and 
the second-moments are predetermined to solve the transport equations for 


u^UjU^,. That Is the values of the mean velocities, u^u^, k, and c 
computed with the method described in the preceding chapter are stored, and 
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where 


+ a 4 (a ik 8 jt,t + a 1J a ki,l + a jk a 1l,t )} 


a 1j “ Vj ‘ 3 4 1j k 




ij,k ax^ 


( 20 ) 


The coefficients are given In Table 2. 

TABLE 2. Coefficients for a\ 


•1 

°2 

®3 

a 4 

-8.14x10-3 

-1.72x10-2 

-4.80x10-2 

-1.02X10- 1 


■i 


.X 



these values are retrieved when the transport equations for u^UjU^ are 
coreputed. In this way, the existing models for the triple-velocity products 
are also evaluated and compared with the present third-moment closure model. 

Since Eq. (15) represents the transport equations of uuu, uuv, uvv, and 
vvv in a two-dimensional coordinate system, and all four components are asso- 
ciated among themselves In their own equations, these four equations are 

solved Iteratively. The Iteration was terminated when the relative residual 

-12 

source of each equation dropped below 3x10 

Figures 12-14 compare the triple-velocity product distributions for dif- 
ferent values of the coefficient C p : 5.0, 6.0, and 10.0. The results are 

5 

also compared with the experimental data of Chandrsuda and Bradshaw. As Is 

observed in this figure, the results are relatively sensitive to the value of 

C p . As Cp becomes larger the levels of the triple-velocity products 

decrease since this term (Eq. (14)) acts as a sink of u^u^u^. After a 

number of parametric tests It Is found that the value of Cp * 5.8 gives the 

23 

best agreement with both the experimental data of Driver and Seegmlller 

5 

and of Chandrsuda and Bradshaw. This optimum value corresponds ;o an 
Inverse coefficient C p ■ 1/Cp * 0.17 which Is slightly larger than the 

I 

value of Cp - 0.11 recommended for the Hanjalic and Launder model (see 
Eq. (17)). 

Figures 15 - 20 compare the present computations of the triple-velocity 
product profile with the experimental data of Chandrsuda and Bradshaw and of 
Driver and Seegmlller, along with the four algebraic models at several dif- 
ferent locations downstream of the step. It is observed that two peaks appear 
across the shear layer: one along the core of the separated shear laypr, the 

other near the bottom wall. Although all the computations exhibit a similar 
trend to the experimental data, the predicted levels of u^UjU^ arc* 


f 


fairly different. While the node! of Shir gives the lowest levels for all 
three components, the model of Cornack et al. provides high levels for uvv and 
vvv but not for uuv. Unlike the model of Cormack et al., the Daly-Harlow 
model gives high levels only for uuv and relatively low levels for both uvv 
and vvv. Among the algebraic models the results with the Hanjallc-Launder 
model seem to be consistent showing nearly the same levels for all three com- 
poiients. This Is because the Hanjallc-Launder model is the only one that has 
a symmetric property In all three directions among the algebraic correlations. 

Agreement between the present computations with the transport equations 
and the experimental data Is generally much better than that using the alge- 
braic models for two reasons. First, the coefficient, C p , can be suitably 
adjusted In the transport equations whereas the algebraic models have been 
tested only for relatively simple shear flows and r.ot for the reattaching 
shear layer. Second, the transport equations of u^u^ have a symmetric 
property in all three directions; thus, the prediction with the transport 
model is better than the prediction with algebraic models for Inhomogeneous 
flows as well as for homogeneous ones. However, the most Important point we 
should note from these results is that the location of the peak along the 
separated shear layer Is predicted more accurately by using the transport 
equations than by using algebraic models. This Is primarily because the 
transport model takes both the convection and the generation due to mean 
strains into account; thus, the trail of energy propagation for the 
triple-velocity products accord'; with the experimental data. 

The results in Figs. 15-2C also suggest that the transport model needs to 
be Improved to give better profiles near the wall. It Is shown in these 
figures that the computed levels are much higher than the experimental data 
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near the wall. In the wall adjacent region the model must Incorporate a low- 
Reynolds number effect although the computations In the free shear flow region 
agree well with the data. The effect of the wall boundary condition on the 
profiles of UjUjU^ near the wall was tested by adopting several other 
models for the numerical nodes next to the wall. However, the Influence of 
the wall boundary condition was negligibly small. Therefore, to improve the 
prediction of the triple-velocity products In the near-wall region. It is 
advisable to develop a low-Reynolds number model In order to take the viscous 
effect into account. 

3.4 Summarizing Remarks 

The development of the third-moment closure has been pursued through this 
project. It is observed that the behavior of the triple-velocity products Is 
completely different from that In the boundary layer flows showing rapid decay 
toward the solid wall in the separated shear layer. 

It Is noted that the transport equation model dramatically Improves the 
prediction of the triple-velocity products throughout the whole flow region 
when compared with the algebraic third-moment models. However, more extensive 
tests are required for Improving the transport equation model since the clos- 
ing of the diffusion and pressure-stress terms of the transport equation has 
not yet been justified. Moreover, it is recomr, ended that a low-Reynolds 
number near-wall model be developed for the transport equation model. 



4. THE ALGEBRAIC-STRESS MODEL 


4.1 Preliminary Remarks 

While It Is Important to Improve the prediction of the Reynolds stresses 
in a wake region. It Is equally Important to develop a better expression for 
empirical coefficients for the diffusion rate of the momentum equation. This 
Is because empirically determined coefficients would Increase uncertainties 
and, thus the prediction becomes less accurate for more geometrically complex 
turbulent flows. 

The preceding chapters have dealt with the Reynolds-stress model. How- 
ever, the algebraic-stress and the k - c models, which are the simpler 
models of turbulence, are considered In this chapter In order to investigate 
the formulation of the coefficient by separating from many unknowns 
which exist In the Reynolds-stress model. If the Reynolds-stress model (RSM) 
were also used, it would be rather difficult to differentiate between the 
errors due to the RSH and the formulation. For futher simplification 
of the problem, the solid wall is completely eliminated in the wake region, 
thus considering only the flow behind a disk. 

The development of Is based on the algebraic-stress model (ASH) 

g 

which was originated by Rodi. Although the Reynolds-stress transport model 
as described in the preceeding chapters Is one of the most advanced model, the 
algebraic-stress model (ASH) Improves the prediction over the k - e model. 
Probably the most advantageous point of the algebraic-stress model (ASH) is 
Its simplicity; the solution procedure of a transport equation is not re- 
quired. The Reynolds stresses are determined by merely solving algebraic 
equations without using the line relaxation or tr^dlagonal matrix solvers, and 



thus It Is much more economical In comparison with the Reynolds-stress 
transport model. 

The above mentioned algebraic-stress model (ASH) Is used to formulate the 

diffusion coefficient of the momentum equations. The model Is compared with 

the conventional empirical constant for the computations of wakes. Since the 

formulation has to be general, axlsymmetric flows have been chosen so that 

2 2 2 

three normal components of the stresses (u , v , and w ) and the shear 
stress (uv) are solved In the region behind a disk. 

4.2 Formulation 

The ASH for the axlsymmetrlc flow can be obtained by using the propo:al 
of Rodi. Rodi proposed that the convection-diffusion of the k-equation Is 
proportional to that of the Reynolds-stress transport equation such as 


U 1 U 1 ak au 1 u i 

k (' U ml^- D < k » “' U ml^- D < u i u j> 


( 21 ) 


where D(k) is given as 


a r ak „ k ak , 

D(k) “ ssg [v 55J + c k ' c u i u j ^ 


( 22 ) 


and where D(u^Uj) represents the diffusion rate of u.jUj. The final 
form obtained is given as follows: 


U 1 U 1 2 


3 D 1j 


H-W, -1*1,6) ♦C W 2^(6 11 - D^) 

r k 3/2 

' C 1 - C wl V* 


(23) 


The ASH is comprised of four algebraic equations for the four Reynolds 
stresses for a two-dimensional axisymmetric flow. These arc given as follows: 
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(26) 


uv « k 


3/2 

(1 - C ? ) G 1? •»• C w? ~ (G, ? - D u ) 
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[G . q _ r L_ 
[ C + C 1 ‘'Hi cy 
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The genera 


tion rate of turbulence kinetic energy (Eq. (7)) Is defined as 


G * -[uv ( |U + |X )+u 2|M +v 2|M + u 2y3 


(28) 


.2 ..2 


The components of the generation (Eq. (7)) for Reynolds stresses (u , v , 
w 2 and uv) are given as 




Gn --2 (u 2 ^ + uv f r ) 


(29) 


2 3V — 3V. 

C 22 ■- 2(v 57* “ v 57> 


(30) 


2 V 

G 33 “ * 2 * r 
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G 12 .-2(v 2 ^ + u 2 |J-uv jf) 


(32) 


The pressure-strain correlation as given by Eq. (6) can be written as: 


♦ll “ " C 1 k 


C, h <u 2 - | k) - C 2 (G n - | G) + [C wl f (u 2 - | k) 


C w2* 2uv 
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The secondary generation rates (Eq. (7)) are given as: 
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Tho model described above Is also compared with the k - e tsodcl. In the 
k - c model the Reynolds stresses are confuted by using the Bousslnesq 
viscosity concept. That Is. all the Reynolds stresses are expressed as 
follows: 


.2 . I 


. A uu 

pk - 2 57 


2 2 . _ av 

pv - 3 pk - 2 v t 3r 


2 2 v , V 

pw ** j pk ~ 2 r 


» sjj av x 
v t ( ar ax' 


Here v t Is evaluated by using the following equation: 


v t « pk /e 


In this study four different models are used out of which two are the 
k - c models and the other two are ASMS. The difference In the two k - c 
models depends upon the value given to C^, whether It Is a constant value 
or a functional value derived from the ASM formulation. The difference in the 
two ASHs Is also the same. Hence, we can define the four models as follows: 


ASH I 


stresses computed by ASM: ■ 0.09 


ASM II — stresses computed by ASM: - f 

k-c I (standard k - c model) -- stresses computed by k - c: - 0.09 

k-c II (modified k-c model) -- stresses computed by k - c: C » f 



» 


The function f is derived upon taking the recirculating effects due to 
two-directional mean strains into account in the formulation. The final form 
is given as 


1 - C 
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The details of the formulation are given in Appendix. 

The constants used In the equations above ere given In Table 3. 

TABLE 3. Constants 


c k 

c wl 

c w2 

Cl 

c 2 

0.1 

0.125 

0.015 

1.5 

0.4 

4.3 Computed Results and 

Discussion 





4.3.1 Computing Details 

Exploratory tests were made for different mesh sizes to Investigate an 
optimum grid Independent state. Figure 22 gives mean velocity and the 
Reynolds-stress profiles for different grid sizes at x/D ■ 6 for the flow 
geometry In Fig. 21(a). 

It Is clear that the velocity and the Reynolds-stress profiles approach 
an equilibrium state for a grid density higher than 42x42. The percentage 
increase In the properties from one grid size to another Is given In Table 4. 
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TABLE 4. Grid test at r/D = 0.5 


Hesh Sizes 

32x32 to 42x42 

42x42 to 52x52 

52x52 to 62x62 

U 

55. OX 

40. OX 

5. OX 

U 2 

75. OX 

25. OX 

o.ox 

v2 

66. 6X 

33. 4X 

O.OX 

w2 

71. 5X 

28. 5S 

O.CX 

uv 

62.5% 

37. 5X 

O.OX 


In Tabic 4 It Is observed that there Is no appreciable change In stresses 
from 52x52 to 62x62, and this Is the reason why we opted for 62x62 as the 
optimum grid size for all the computations In this chapter. 

The Initial values for the Reynolds stresses are given by the solutions 
of the k - e model to ensure stability In the course of the Iteration pro- 
cess for the elliptic equations. The computation Is -.erminated when the 
maximum residual source of the transport equation falls below 0.3 percent of 
the total source of Its equation. Although iteration number depends on many 
factors such as grid size, flow geometry, boundary conditions. Initial 
conditions, Reynolds number, etc., the typical number ranges from 500 - 800 
for the grid size 62x62. 

The computational domain for the external flow past a disk consists of 30 
step height lengths In the x-dlrec'lon and 5 step heights In the r-dlrectlon. 

A computation was also made by extending the length from 30 to 40 step 
heights. When the results obtained with two different geometries were com- 
pared with each other, no appreciable difference was observed in any proper- 
ties. Thus, the geometry of 30 step height lengths is used in x-dlrection for 


all the computations of the external flows. The expansion ratio of 1.05 In 
both x- and redirections Is found to be the best for this flow geometry. The 
computational domain of the second geometry shown In Fig. 21(b) consists of a 
length of 24 step heights In x-dlrcctlon and two step heights In r-dlrectlon. 
The optimum expansion ratios for this case are found to be 1.05 In x and 1.00 
In r-dlrectlons. 

4.3.2 External Flows 

Figure 23 represents the mean velocity profiles at several locations In 

the flow field beyond the body for the flew geometry which Is shown in Fig. 

21(a). These are computed by employing k-e I, k-c II, ASH I and ASM II. 

These results are also compared with the experimental data obtained by 
21 

Carmody. The velocity profiles obtained with the models mentioned above 
agree with the experimental data outside the separated shear flow region but 
some disagreement Is observed In the recirculating region which might be 
attributed to the sharp edge of the disk used in the experimental setup caus- 
ing a larger recirculation region than that predicted in the present study. 

Figures 24-27 represent the Reynolds stresses obtained by the ASM and the 
k-c model which are compared with the experimental data. The agreement of 
the stresses obtained by ASH I with the experimental data seems to be better 
than agreement obtained by the other models. The other models have a similar 
trend but do not agree with the experiment as closely as ASH I. The agreement 
o'; all the numerical models with the experimental data Is better In the region 
outside the recirculation region. Since the k-c models (k-c I and k-c 
II) cannot take the non-lsotroplc effect into account, the errors In the 
results obtained by the k-c models are significant, particularly in the 
recirculating region. Despite the improvement with the use of the ASM models. 
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these models produce a second peak in the recirculating region (near the flow 
Impingement point), which Is caused by the streamwlse acceleration of the flow 
near the flow Impingement point. This Is verified by a test In which the 
gradient of the mean velocity In x-dlrectlon ( ati/ax) Is omitted resulting 
In diminishing of this second peak. 

From the results mentioned above It Is concluded that the algebraic- 
stress models provide better results In the whole domain of the flow geometry 
than the k - e model. This Is evident If we observe the results around the 
flow impingement point (x/D ■ 3.0) as summarized In Table 5. 

TABLE 5. Percentage error for Reynolds stress at r/D « 0.5 In comparison 
with experimental data 




uv 


u 2 

V 

2 

w 2 


ASH I 

k-c I 

ASH I 

k-e I 

ASH I 

k-e I 

ASH I 

k-c I 

x/0-3 

30% 

65* 

51* 

66* 

30* 

36* 

4* 

33* 

x/D-6 

21* 

71* 

23* 

61* 

29* 

50* 

34? 

53* 

x/D-9 

10* 

13* 

10* 

40* 

25* 

54* 

18* 

41* 


Figure 28 shows the turbulence kinetic energy profiles obtained by using 

the models considered In this study which are compared with the experimental 

data. These profiles also show that ASH I Is better than the other models. 

Figure 29 shows the C distribution at different locations In .the flow 

u 

field and Fig. 30 shows the G/c profiles at the same locations. It Is noted 
that the profiles for ASH II and k-c II are very close to each other 
and so are the results computed with them. It Is also found that there Is a 
strong dependence of stress levels on the G/c distribution. If G/c is 
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higher, then Is lower Improving the Reynolds-stress levels; thus ASH I, 

which has the smallest C In most of the flow domain, Is giving the best 

v 

prediction of the Reynolds stresses among the models considered. The position 
of the peak of G/c profiles shown In Fig. 30 computed by ASH I Is higher 
than the peaks produced using the other models. This Is because ASH I creates 
only a small spreading rate. 

4.3.3 Internal Flows 

For the computations of Internal flows only ASH I and k-c I are em- 
ployed instead of using four models. Figure 31 represents the mean velocity 

profiles In the flow direction computed by using ASH I and k-c I. These are 

22 

compared with the experimental data of Taylor and Whltelaw for D/Dy « 

0.5 (Fig. 21(b)). The results obtained by using ASH I agree well with the 

experimental data and are better than those obtained using k-c I. Figures 

2 2 2 

32-34 represent the three Reynolds-stress profiles (u , v and w ) 

obtained by using k-c I and ASH I which are compared with experimental 
22 

data. In these figures it Is also shown that ASH 1 gives better results 
for the flows In this geometry. 

Even In the internal flow geometry the second peak In the Reynolds-stress 

profiles appears in the recirculating region near the flow impingement point 

when using ASH I. This behavior can be explained in the same manner as done 

in the case of the external flow geometry. 

The stresses shown in Figs. 32-34 are very small outside the mixing layer 

because of negligible generation -ates. The stresses then grow slightly near 

the top wall where the mean velocity decreases rapidly due to the viscous 

2 2 

effect from the wall. It is observed that the stresses u and w' are 
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higher than v^; this feature Is similar to that of the wall boundary layer 
flows. 

Figure 35 shows the turbulence kinetic energy distribution obtained with 
both k-e I and ASH I, and they are compared with the experimental data. The 
turbulence kinetic energy shows the highest value near the flow Impingement 
point and It decays towards downstream. Here also it Is shown that the re- 
sults computed by using ASH I agree better with the experimental data. Figure 
36 shows the G/c distribution obtained by ASH I and k-c I. The trend Is 
similar to the G/c distribution for the external flow (see Fig. 30). 

4.4 Summarizing Remarks 

The Investigation In this chapter revealed that even the algebraic stress 
model Improves the prediction of the Reynolds stresses considerably when com- 
pared with the k-e model. However, It Is also observed that the formula- 
tion of based on the algebraic stress model does not Improve the 
prediction at all. This observation suggests that any model based on the 
Bousslnesq viscosity concept never Improves the computations of the Reynolds 
stresses. For this reason, we should discard the k-e model type of 
turbulence model and should create second- or higher-moment closure models for 
the computation of recirculating flows. 
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5. APPLICATION TO HEAT TRANSFER PROBLEMS 

5.1 Preliminary Remarks 

The models developed In the course of this project can be applied to cany 
Industrial problems without making many modifications. One application we 
exercised recently was the computation of the flow and heat transfer rates In 
periodically corrugated wall channels which appear In a corrugated-wall- 
channel heat exchanger. The Reynolds-stress model was primarily employed for 
the computations of the hydrodynamic turbulence variations and the heat- 
momentum analogy was used to evaluate the heat transfer rates. The results 
were compared with the computations made by using the k - c model. The main 
concern regarding the differences produced by the two different closure models 
is whether the non-lsotroplc effect in the corrugated wall channel Is predomi- 
nant or not. If It Is, a significant difference will occur In the results of 
the computations. 

5.2 Computational Method 

The computational domain of the flow field Is the region ArCOEJIHGF (in 
Fig. 37). This region consists of two cycles (ABCOIHGF) and the additional 
region (OEJI). For computations of laminar flows, only one cycle of the 
corrugated wall channel may be used with a periodic inlet/outlet condition; 
however, a two-cycle system Is needed for complex turbulent flow computations 
because of numerical instability. In turbulent flow equations more nonlinear 
transport coefficients appear which results In a very poor hlsto.-y of conver- 
gence performance in the course of the Iteration process. In addition, a 
two-pass procedure is employed. That is, the computation was started out with 
a constant prescribed Inlet condition at AF and an outflow condition at EJ. 
With a line relaxation method, computations were iterated until relative 
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residual sources of all the transport equations decreased below 28. Then the 
periodic condition was activated by transferring the updated values at the 
section 01 to the Inlet section AF. Finally, the computation was terminated 
when all the relative residual sources became less than 0.13. 

At the wall boundaries, the momentum, energy and the turbulence kinetic 
energy wtrc evaluated from the "law of the wall" while the energy dissipation 
rate was determined from the "local equilibrium condition," Which gives linear 
variation of the turbulence length scale from the solid wall. 

The average friction factor Is computed as 



where L is the distance between corresponding planes and aP is the average 
pressure drop between these planes. 

The local Musselt number is obtained by the following equation 


Nu *= 


(q w /r)2b 

' T w - V 


(48) 


where the constant heat flux condition Is employed. The bulk temperature, 
T^, is determined as 


T . mui-dy 

b J |U|dy 


( 49 ) 


where the Integrals are to be carried over the cross-sectional area of the 
channel. The absolute value of thr; velocity is taken so that the regions with 
reverse flows are also properly represented. Accordingly, the averaged 
Husselt number is defined as 
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where x Is the coordinate along the channel wall. 

5.3 Computed Results and Discussion 

5.3.1 Validation Test of the Reynolds Stresses 

Due to lack of experimental data for corrugated wall channel flows, the 
present numerical model cannot be justified for all the variables obtained. 

In particular, turbulence quantities such as the Reynolds stresses need to be 
compared with some experimental data. For this reason, data for a backward 
facing step flow was chosen to perform a validation test of the Reynolds 
stresses since the wall region behind a step Is similar to the flow field 
along the wall ABCDE In Fig. 37. 

Figure 38 shows the computed results of the mean velocity and the 

Reynolds-stress distributions In the region behind the step. These results 
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are compared with the experimental data of Driver and Seegmlller. Here 52 
x 52 grid points are used for the computational region of 50a x 5a. Both in 
recirculating and recovering regions, agreement between the present computa- 
tions and the experimental data Is reasonably good (within 30%). 

5.3.2 Average Nusselt Number and Skin Friction Coefficients 

0 3 

Figure 39 shows Nu/Pr ’ as a function of Reynolds number. The 

24 ?5 

experimental data of both Izumi et al. and of O'Brien and Sparrow* are 
used for comparison with the present computations. Although the geometry 
considered by O'Brien and Sparrow has a bend of 120° and a/b * 1.6, the levels 
of the average Nusselt number seem to be almost the same as those given by 
Izumi et al. In which a channel with a 90° bend and a/b <= 2 were used. 
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For a/b - 2, the confutations by tha k - e model and the present 
model (RSH) are also compared. It Is noteworthy that the RSH Improves the 
prediction by 30% compared with the k - e model. This observation Is 
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consistent with the case of an Infinitely long channel with two bends In 
which It was also discerned that the prediction of the Nusselt number was 
Improved by 20-30% by employing the RSH. This Is because the turbulence 
levels In recirculating regions are correctly evaluated since the RSH accounts 
for non-isotropic behavior of the Reynolds stresses. 

In Fig. 39 the computations for a/b - 1.5 and 1.25 are also shown. It Is 
observed that the difference In the Nusselt number created by the change In 
corrugation period Is minor showing only a 2-3% decrease from a/b = 1.5 to 
2.0. This result is also consistent with the experiments of both Izumi et al. 
(a/b =2) and O'Brien and Sparrow (a/b « 1.6). 

Figure 40 shows the skin friction coefficient as a function of Reynolds 
number. Unlike the case of the Nusselt number, the skin friction rate depends 
on the corrugation period, a/b. It Is observed that the dependency of f on Re 
is constant for both a/b = 1.25 and 1.5, but for a/b «= 2 it decreases slightly 
as the Reynolds number increases for Re > 3000. In Fig. 40 experimental data 
for a straight channel is compared with the computations. This figure shows 
that the computed curves approach the straight channel data as a/b decreases 
to zero. 

5.3.3 Velocity and the Reynolds-Stress Profiles 

Figure 41 demonstrates the velocity vectors in different cl annals. As is 
shown, recirculating flow is created at each concave corner of the channel. 

The flow Impingement position and the recirculating region are dependent on 
the corrugation period, a/b. 


Computed mean velocities and the Reynolds stresses In the channel where 
o/b - 2 and for Re - 3000 are shown In Fig. 42. At section B, the mean veloc- 
ity becomes maximum near the wall AB, and It decreases sharply toward the wall 

2 

AB. This high shear stress Is causing high levels of the stresses u and 
2 

v In that region. If we look at the region near the opposite side of the 

channel (FG wall), the mean velocity profile Is rather smooth although the 

velocity becomes negative near the wall FG due to a recirculating flow; thus, 

the normal stresses are not as high as the shear stress uv. The levels of 
2 2 

u and v Interchange In the next bend section (cross-section CH) showing 

2 2 

a relatively high level of v and a low level of u because the y com- 
ponent of velocity Is predominant over the x component mean velocity. At 

2 2 
section H, the u level is higher near the wall HI, whereas the v level 

Is high near the opposite side of the channel where the y component of the 

velocity Is still large due to the flow deflection caused by the separation at 

2 

the corner H. The maximum point of u moves to the wall DE as we move 
toward downstream. 

5.3.4 Flow in a Channel with Fins 

Computations were made for flows in channels which have fins inserted at 
each bend to cause more deflection of the flow. Local Nusselt number distri- 
butions are shown in Figs. 43 and 44 for a/b = 1.25 and 1.5, respectively. 
Although the general pattern of the distribution is similar for both a/b « 

1.25 and 1.5, the distribution depends on the length of the fin. When fins 
are Inserted, the position of the peak of the Nusselt number moves toward the 
corner C along the wall BC. This Is mainly due to the expansion of the 
recirculating region downstream of the corner B with the insertion of the 
fins. A similar trend also appears near the wall D. It is commonly observed 




that the heat transfer rates Increase as the length of the fin Increases up to 
t fifj /b - 0.5. However, the case of L f1fJ /b ■ 0.75 does not display propor- 
tional Increases In the heat transfer rates. 

The dependence on the length of fins Is also shown In Fig. 45 for the 
average Nusselt number. Here the average Nusselt number Increases as L f1n 
Increases from 0 to 0.5b, but It stays at constant values for longer fins. 

The rate of Increase In the Nusselt number Is larger for higher Reynolds num- 
bers. 

5.4 Summarizing Remarks 

In this chapter, we saw an application of the Reynolds-stress model to a 
heat exchanger problem. One of the most Impressive features shown here Is 
that the Reynolds-stress model improves the computation of the heat transfer 
rates by 20-30X. This Is mainly because the heat transfer rates are almost 
proportional to the hydrodynamic variables, thus If the turbulence levels are 
evaluated correctly so are the heat transfer rates. 



6. FINAL REMARKS AND FUTURE V50RK 




In this project both the second- and third-moment closures have been In- 
vestigated. The models developed revealed significant Improvement In the 
prediction of complex turbulent shear flows that are accompanied by reattach- 
ment, recirculation, and flow deflection. However, the tests were made only 
for each component of each model— that Is, the third-moment transport equa- 
tions were solved without coupling to the second-moment equations. Although 
each modeling seems to be successfully carried out, there Is no guarantee that 
these models perform equally well when all the transport equations are solved 
simultaneously. At that stage It Is anticipated that a huge amount of 
computational time and storage will be required since numbers of transport 
equations have to be solved by employing sufficiently high density of mesh 
size. Probably the problem. If any, will be numerical Instability rather than 
a difficulty due to mathematical complexity. For this reason we believe it is 
extremely Important to complete the component study pt ior to consolidating all 
the models in one package of the computer code. 

It should also be noted that there is plenty of room for Improvement of 
the presently developed second- and third-moment models. For example, a low 
Reynolds number model coupled with a near-wall correction must be incorporated 
for the third-moment closure equations. Secondly, the development of the 
equivalent models for turbulence scalar quantities (heat or chemical species) 
will undoubtedly advance the state-of-the-art technology in turbulence model- 
ing for applications of many industrial problems. 



APPENDIX 


OF A$fl f Algebraic Stress Mnrf^) 

The transport equation for the kinetic energy can be written as 
Ot - D(k) + G - c 


The transport equation for the Reynolds stresses Is given by 


Dt ‘ D <Vj) + G tJ -u + ^ (A . 

where D(u lUj ) Is the diffusion rate of the Reynolds stresses. Rodl 8 

proposed that the convection-diffusion string of the k-equatlon Is assumed to 

be proportional to that of the Reynolds-stress transport equation, as given In 
the equation below. 
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From (A. 2) and (A. 3) we have 
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Now, substituting Eq. (A.l) in (A. 4) we have 
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where Is defined In Eqs. (33)-(36). 
Daly-Harlow, 17 we have 
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Using the above result In (A. 5) we have 
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After rearranging, we have 
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Regrouping the terns we have 
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Finally, the equation for Reynolds stresses Is given as follows: 
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By neglecting wall effects and the secondary mean strains, u 2 , v 2 and uv 
are given as follows 
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Substituting Eqs. (A.9) and (A.10) Into (A.ll) and simplifying. we get 
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After multiplying both sides by p and rearranging, we have 
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Figure 10 Dissipation of turbulence energy 
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Triple-velocity products 






Figure 16 Triple-velocity products 
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Figure 18 Triple-velocity products 


HI o 









x/H=6.0 9 x/H=8.0 O x/ 















o EXPERIMENT (Cirmody) 

ASM I 

ASM a 

k-£ I 

k£fl 



U 2 /U2 


figure 24 


Streamwise component of Reynolds -stress distributions for the 
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Figure 26 Azimuth component of Reynolds-stress distributions for the ex 
ternal flow over a disk 
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Figure 41 Velocity vectors in channels 
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